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Abstract. We have used high-resolution, HST WFC3/IR, near-infrared imaging to conduct a 
detailed bulge-disk decomposition of the morphologies of ~ 200 of the most massive (M* > 
ii . 10 11 M ) galaxies at 1 < z < 3 in the CANDELS-UDS field. We find that, while such massive 
galaxies at low redshift are generally bulge-dominated, at redshifts 1 < z < 2 they are predom- 
inantly mixed bulge+disk systems, and by z > 2 they are mostly disk-dominated. Interestingly, 
^ ■ we find that while most of the quiescent galaxies are bulge-dominated, a significant fraction 
(25 — 40%) of the most quiescent galaxies, have disk-dominated morphologies. Thus, our re- 
suits suggest that the physical mechanisms which quench star-formation activity are not simply 
connected to those responsible for the morphological transformation of massive galaxies. 
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. 1. Introduction 

The study of the high- redshift progenitors of today's massive galaxies can provide us 
with invaluable insights into the key mechanisms that shape the evolution of galaxies in 
the high-mass regime. 

The latest generation of galaxy formation models are now able to explain the number 
densities and ages of massive galaxies at high redshift. However, this is only part of the 
challenge, as recent studies have posed new questions about how the morphologies of 
massive galaxies evolve with redshift. 
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In addition to the basic question of how high-redshift galaxies evolve in size, there is 
also still much debate about how these massive galaxies evolve in terms of their fun- 
damental morphological type. Extensive studies of the local Universe have revealed a 
bimodality in the colour-morphology plane, with spheroidal galaxies typically inhabiting 
the red sequence and disk galaxies making up the blue cloud (e.g. Baldry et al. 2004). 
However, recent studies at both low (e.g. Bamford et al. 2009) and high rcdshift (e.g. van 
der Wei et al. 2011) have uncovered a significant population of passive disk-dominated 
galaxies, providing evidence that the physical processes which quench star-formation may 
be distinct from those responsible for driving morphological transformations. This result 
is particularly interesting in light of the latest morphological studies of high-redshift 
massive galaxies by Buitrago et al. (2013) and van der Wei et al. (2011) who find that, 
in contrast to the local population of massive galaxies (which is dominated by bulge 
morphologies), by z ~ 2 massive galaxies are predominantly disk-dominated systems. In 
this work we attempt to provide significantly improved clarity on these issues. 

2. Data and Morphological Fitting 

The CANDELS (Grogin et al. 2011, Koekemoer et al. 2011) near-infrared F160W 
data provides the necessary combination of depth, angular resolution, and area to enable 
the most detailed study to date of the rest-frame optical morphologies of massive (M* > 
10 11 M©) galaxies at 1 < z < 3 in the UKIDSS Ultra Deep Survey (Lawrence et al. 2007). 
For this study we have constructed a sample based on photometric redshifts and stellar 
mass estimates which were determined using the stellar population synthesis models of 
Bruzual & Chariot (2003) assuming a Chabrier initial mass function (see Bruce et al. 2012 
for full details). This provides us with a total mass-complete sample of ~ 200 galaxies. 

We have employed the GALFIT (Peng et al. 2002) morphology fitting code to de- 
termine the morphological properties for all the objects in our sample. To conduct the 
double component fitting we define three components: a Sersic index fixed at n = 4 bulge, 
an n = 1 fixed disk and a centrally concentrated PSF component to account for any AGN 
or nuclear starbursts within our galaxies. These three components are combined to gener- 
ate six alternative multiple component model fits, of varying complexity, for every object 
in the sample. These models are formally nested, and thus x 2 statistics can be used to 
determine the "best" model given the appropriate number of model parameters. 

3. Evolution of Morphological Fractions 

Armed with this unparalleled morphological information on massive galaxies at high 
redshift we can consider how the relative number density of galaxies of different morpho- 
logical type changes during the key epoch in cosmic history probed here. 

In Fig. 1 we illustrate this by binning our sample into four rcdshift bins of width 
Az = 0.5, and consider three alternative cuts in morphological classification as measured 
by B/T from our bulge-disk decompositions. In the left-hand panel of Fig. 1 we have 
simply split the sample into two categories: bulge-dominated (B/T > 0.5) and disk- 
dominated (B/T < 0.5). In the central panel we have separated the sample into three 
categories, with any object for which 0.3 < B/T < 0.7 classed as "Intermediate". Finally, 
in the right-hand panel we have expanded this Intermediate category to encompass all 
objects for which 0.1 < B/T < 0.9. 

From these panels it can be seen that z ~ 2 marks a key transition phase, above 
which massive galaxies are predominantly disk-dominated systems and below which they 
become increasingly mixed bulge+disk systems. We also note that at the lowest redshifts 
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Figure 1. The redshift evolution of the morphological fractions in our galaxy sample, after 
binning into redshift bins of width Az = 0.5 and using three alternative cuts in morphological 
classification (both to try to provide a complete picture, and to facilitate comparison with 
different categorisations in the literature). 



probed by this study (z ~ 1) it is seen that, while bulge-dominated objects are on the 
rise, pure-bulge galaxies (i.e. objects comparable to present-day giant ellipticals) have yet 
to emerge in significant numbers, with > 90% of these high-mass galaxies still retaining a 
significant disk component. This is compared with 64% of the local M* > 10 11 M galaxy 
population, which would be classified as pure-bulges from our definition (B/T > 0.9, 
corresponding to n > 3.5 ) from the sample of Buitrago et al. (2013). Thus, our results 
further challenge theoretical models of galaxy formation to account for the relatively rapid 
demise of massive star-forming disks, but the relatively gradual emergence of genuinely 
bulge-dominated morphologies. 



4. Star-forming and Passive Disks. 

In addition to our morphological decompositions we also make use of the SED fitting al- 
ready employed in the sample selection to explore the relationship between star-formation 
activity and morphological type. 

Fig. 2 shows specific star-formation rate (sSFR) versus morphological type for the 
massive galaxies in our sample, where morphology is quantified by single Sersic index in 
the left-hand panel, and by bulge-to-total Hieo-band flux ratio (B/T) in the right-hand 
panel. The values of sSFR plotted are derived from the original optical-infrared SED fits 
employed in the sample selection, and include correction for dust extinction as assessed 
from the best fitting value of Ay derived during the SED fitting. As a check of the po- 
tential failure of this approach to correctly identify reddened dusty star- forming galaxies, 
we have also searched for 24 /im counterparts in the Spitzer SpUDS MIPS imaging of the 
UDS, and have highlighted in blue stars those objects which yielded a MIPS counterpart 
within a search radius of < 2 arcsec. 

To first order, our results show that the well-documented bimodality in the colour- 
morphology plane seen at low redshift, where spheroidal galaxies inhabit the red sequence, 
while disk galaxies occupy the blue cloud is at least partly already in place by z ~ 2. 

Nonetheless, the sample also undoubtedly contains star- forming bulge-dominated galax- 
ies and, perhaps more interestingly, a significant population of apparently quiescent disk- 
dominated objects. To highlight and quantify this population we have indicated by a box 
on both the panels the region occupied by objects with disk-dominated morphologies and 
sSFR < 10 _10 yr —1 . In the left-hand panel, disk-dominated is defined as n < 2.5, and 
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40 ± 7% of the quiescent galaxies lie within this box (if we exclude the 24 /im detections), 
while in the right-hand panel, disk-dominated is dchned by B/T < 0.5, in which case 
25 ± 6% of the quiescent objects lie within this region. 
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Figure 2. Plots of specific star-formation rate (sSFR) versus morphological type as judged by 
single Sersic index (left-hand panel) and bulge-to-total _ffi6o-band flux ratio (B/T) (right-hand 
panel) . 

The presence of a significant population of passive disks among the massive galaxy 
population at these redshifts indicates that star-formation activity can cease without 
a disk galaxy being turned directly into a disk-free spheroid, as generally previously 
expected if the process that quenches star formation is a major merger. 

One possible mechanism for this arises from the latest generation of hydrodynami- 
cal simulations (e.g. Keres et al. 2005, Dekel et al. 2009a) and analytic theories (e.g. 
Birnboim & Dekel 2003), which suggest a formation scenario whereby at high redshift 
star-formation is fed through inflows of cold gas. Another scenario which can account 
for star-formation quenching, whilst still being consistent with the existence of passive 
disks, is the model of violent disk instabilities (e.g. Dekel et al. 2009b), coupled with 
"morphology quenching" (Martig et al. 2009). 
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